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Abstract. In this paper we develop a symbolic technique to obtain asymptotic expressions 
for ruin probabilities and discounted penalty functions in renewal insurance risk models when the 
premium income depends on the present surplus of the insurance portfolio. The analysis is based on 
boundary problems for linear ordinary differential equations with variable coefficients. The algebraic 
structure of the Green's operators allows us to develop an intuitive way of tackling the asymptotic 
behavior of the solutions, leading to exponential-type expansions and Cramer-type asymptotics. 
Furthermore, we obtain closed-form solutions for more specific cases of premium functions in the 
compound Poisson risk model. 
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1. Introduction. The study of level crossing events is a classical topic of risk 
theory and has turned out to be a fruitful area of applied mathematics, as (depending 
on the model assumptions) often subtle applications of tools from real and complex 
analysis, functional analysis, asymptotic analysis and also algebra are needed (see e.g. 
[4] for a recent survey). 

In classical insurance risk theory, the collective renewal risk model describes the 
amount of surplus U{t) of an insurance portfolio at time t by 



where c represents a constant rate of premium inflow, N{t) is a renewal process that 
counts the number of claims incurred during the time interval (0,t] and {Xk)k>o is a 
sequence of independent and identically distributed (i.i.d.) claim sizes with distribu- 
tion function Fx and density fx (also independent of the claim arrival process N(t)). 
Let (Tfc)fc>o be the i.i.d. sequence of interclaim times. One of the crucial quantities to 
investigate in this context is the probability that at some point in time the surplus in 
the portfolio will not be sufficient to cover the claims, which is called the probability 
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of ruin 

i;{u) = P{Tu<oo\U{0)^u), 

where U{0) = u > is the initial capital in the portfoho and 

Tu = inf {t>0: U{t) < | C/(0) = u}. 

A related, more general quantity is the expected discounted penalty function, which 
penalizes the ruin event for both the deficit at ruin and the surplus before ruin, 

= E (e-*^" w{U{Tu~), \U{Tu)\) 1t„<oo | C/(0) = u) , 

where 5 > is a discount rate and the penalty w{x, y) is a bivariate function. 
is often referred to as the Gerber-Shiu function, see [8]). 

The classical collective risk model is based on the assumption of a constant pre- 
mium rate c. However, it is clear that it will often be more realistic to let premium 
amounts depend on the current surplus level. In this case, the risk process (1.1) is 
replaced by 



U{t)=u + / p{U{s)) ds-Y^Xk 



Hence, in between jumps (claims) the risk process moves deterministically along the 
curve ip{u,t), which satisfies the partial differential equation 

at ou 

There are only a few situations for which exact expressions for are known for 

surplus-dependent premiums. One such case is the Cramer-Lundberg risk model 
(where N{t) is a homogeneous Poisson process with intensity A) and the linear pre- 
mium function p{u) = c + su, which has the interpretation of an interest rate s on 
the available surplus. In the case of exponential claims, it was already shown by [21] 
that the probability of ruin then has the form 

"^^"^ " ^A/e^A/s eM-He) + A£Ve-ir(f , A) r(^^^, f ), (1.2) 

where T{ri,x) = t^~^e~*dt is the incomplete gamma function (for extensions to 
finite-time ruin probabilities, see [11, 12] and [3]). In fact, for the Cramer-Lundberg 
risk model with exponential claims and general monotone premium function p(u), one 
has the explicit expression 

■ip(u) = [ -7--T- exp{Ag(x) — fix} dx, (1-3) 

Ju P{X) 

where I/70 = 1 + A p{x)~^ exp {Xq{x) — ^x} dx and q{x) = dy is assumed 

finite for x > (see [22]). Since for surplus-dependent premiums the probabilistic 
approach based on random equations does not work, and also the usual analytic 
methods lead to difficulties because the equations become too complex, it is a challenge 
to derive explicit solutions beyond the one given above. 
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In this paper we will employ a method based on boundary problems and Green's 
operators to derive closed-form solutions and asymptotic properties of ^j{u) and $(w) 
under more general model assumptions. 

For that purpose we will employ the algebraic operator approach developed in [2] . 
However, since that approach was restricted to linear ordinary differential equations 
(LODEs) with constant coefficients, we will have to extend the theory to tackle the 
variable-coefficients equations that occur in the present context. 

In Section 2 we derive the boundary problem for the Gerber-Shiu function in 
a renewal risk model with claim and interclaim distributions having rational Laplace 
transform. For solving it, we employ a new symbolic method, described in Section 3. 
This allows to construct integral representations for the solution of inhomogeneous 
LODEs with variable coefficients, for given initial values, under a stability condition. 
In Section 4 we derive a general asymptotic expansion for the discounted penalty 
function in the renewal model framework. Subsequently, Section 5 is dedicated to 
the more specific case of compound Poisson risk models with exponential claims, 
for which we have second-order LODEs. More specifically, in 5.1 we derive exact 
solutions for a generic premium function p{u). Further, in 5.2, we consider some 
interesting particular cases of p{u). In 5.3 we identify the necessary conditions a 
premium function should satisfy such that the asymptotic analysis is possible and the 
assumptions necessary for the asymptotic results in Section 4 are validated. We will 
end by giving concrete examples of such premium functions and their asymptotics. 

Throughout the paper we will assume that U{t) — > oo a.s. This assumption is 
satisfied for example when p{u) > EX/Et + for some > and sufficiently large u; 
see e.g. [4]. 

2. Deriving the boundary problem. Assume that the distribution of the 
interclaim time of the renewal process N(t) has rational Laplace transform. For 
simplicity of notation, we assume further that the rational Laplace transform has a 
constant numerator. Then its density fr satisfies a LODE with constant coefficients 

£rij^)Mt)^0 (2.1) 

and homogeneous initial conditions fi''\t) = (fc = 0, . . . , ri — 2), where 

Using the method of [6], we can then derive an integro-differential equation for $(u) 
C; (pi^)-^ - = «o ^ - y) dFx{y) + uj{u)^ , (2.2) 

where £* is the adjoint operator of Cr defined through 

Cl{x) = Cr{-x) = i-x)" + a„_i(-x)"-i + ■ • ■ + Q'o. 

Assume now that the claim size distribution also has a rational Laplace transform, 
so that its density fx satisfies another such LODE 

CxAfx{y)^0 (2.3) 
ay 

with initial conditions fi^\x) = (fc = 0, . . . , m — 2), where 
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Then the integro-differential equation (2.2) becomes a LODE with variable coefHcients 
of order m + n, namely 



T^u) = g{u) (2.4) 



with differential operator 



and right-hand side 



T = Cx[j^\C*Ap{u)^-5]~a^p^ (2.5) 



g{u) = ao Cx{-^)u:{u), 
du 



where uj{u) = w{u^ U ~ u) fx{y) dy. For (5 = and w = 1, Equation (2.4) reduces 
to the well-known equation for the probability of ruin. 

The equations hold for sufficiently regular functions p. In the special case p{u) = c 
one recovers the LODE with constant coefficients whose characteristic polynomial is 
of degree n-\~m and corresponds to Lundberg's equation. It is known that, for <5 > 0, 
this polynomial has m solutions Ui, with negative real part, and n solutions pi, with 
positive real part; see for example [14] and [13]. In [2], we have derived 

$(u) = 71 e^^" + • • ■ + 7,ne'^-" + G.g(«), (2.6) 

where the 7^ are determined by the initial conditions and 

fl- nU nCX) poo 

Gg{u) ^J2J2^^A ^'''^""^^ + / ^'^^""^^ " / ^'"^^1 9(0 (2.7) 
i=i j=i "'0 "'0 

defines the Green's operator for the inhomogeneous LODE (2.4) with homogeneous 
boundary conditions, where 

m n 
k=l,k=ii k=l,k=ij 

The boundary conditions for (2.4) consist of the initial conditions (0)(fc = 0, . . . , m- 
1), determined from the integro-differential equation, and the stability condition 
$(oo) = 0, provided by the model assumptions. 

In analogy to the constant coefficients case, we assume the existence of a fun- 
damental system for equation (2.4) with m stable solutions Si{u) and n unstable 
solutions rj{u). Here a solution f{u) is called stable if f{u) and unstable if 
f{u) — > oo as w — )■ oo. We write ti, . . . ,tm+n for the complete sequence of solu- 
tions si , . . . , s„i , ri , . . . , r„ , and we assume furthermore that the successive Wronskians 
Wk = M^[ii, • ■ • , ifc], for k = 1, . . . ,m + n are all nonzero on the half-line M+ = [0, oo). 
Under these assumptions, the algebraic operator approach developed for the constant 
coefficients case [2] will be extended to the surplus-dependent premium case in Sec- 
tion 3, and the general solution of (2.4) then has the form 

$(u) = 7isi(u) H \-j,nS,n{u) + Gg{u), 



where the 7^ are determined by the initial values and Gg{u) is again the Green's 
operator for the inhomogeneous LODE (2.4) with homogeneous boundary conditions. 
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but this time with non-constant p{u). As a consequence, the representation (2.7) 
is no longer valid, and we will derive a new explicit expression that generalizes it 
(Theorem 3.4). 

Let us complete this section with a remark about how to cheek that the funda- 
mental system has stable and unstable solutions. Roughly speaking, this amounts to 
an asymptotic analysis of the solutions of the homogeneous equation. According to [7, 
Ch.5], one can identify conditions on p{u) that guarantee the existence of such a fun- 
damental system. These conditions specify the structure of the coefficients, namely: 
either they converge (sufficiently fast) to constants — in this case one speaks of almost 
constant coefficients — or they diverge to infinity. The canonical form of (2.4) indi- 
cates of course that the former case applies for our setting here. However, the speed 
of convergence of the coefficients depends crucially on the premium function p{u) . For 
instance, we will show in Example 5.4 that for p{u) = ce^^"", the LODE with almost 
constant coefficients converges to the LODE with constant coefficients given in [2]. 

3. Green's operator approach. In the previous section we have seen that the 
core task for computing the Gerber-Shiu function $(u) is to determine the Green's op- 
erator G for the inhomogeneous LODE (2.4) with homogeneous boundary conditions 
consisting of the initial conditions $^'^'^(0) = (fc = 0, . . . , m — 1) and the stability 
condition $(cx)) = 0. In this section we will present a symbolic method that allows to 
construct G for a generic LODE with variable coefficients and homogeneous boundary 
conditions. In other words, we consider boundary problems of the general type 

r$(u) = g(u), 

(3 1) 

$(0) $'(0) = ••■ = $(™-i)(0) = and $(oo) = 0, 

where T = 1?'"+" + c™+„_i(w) h h ci(u) D + co(u) is a linear differential 

operator with variable coefficients (and leading coefficient normalized to unity) and 
D = Under the conditions described in Section 2 the solution of (3.1) is unique 
and depends linearly on the so-called forcing function g{u). Therefore the assignment 
g h^- $ is a linear operator: the Green's operator G of (3.1). The following fact follows 
immediately from the theory of differential equations. 

Theorem 3.1. The Gerber-Shiu function equals 

$(u) = 7iSi(m) + h 7mSm(w) + Gg{u), (3.2) 

where G is the Green's operator for equation (3.1), and the constants can be iden- 
tified from the initial conditions. 

For describing our new method of constructing an explicit representation of G, 
let us recall how this was achieved in [2] for the special case of constant coefheients 
Ci{u) = Ci. We will use the same notation as there, in particular the basic operators 
A = Jp , B — and the definite integral F = A -\- B = j'^ . Employing the basic 
operators, the crucial idea was to factor the Green's operator as 

G= (-1)"A,, ...A,„Bp, ...Bp„, (3.3) 

where the factor operators are defined by Aa- = e°"-^Ae~°'^ and Bp = e^^Be^P'-^ with 
(Ti and pj as described before. So the strategy was to decompose the problem and 
tackle the stable exponents with the basic operator A, the unstable ones with B. 



6 ALBRECHER, CONSTANTINESCU, PALMOWSKL REGENSBURGER, ROSENKRANZ 



This idea can be carried over to the general case of (3.1). Using the results of 
[19], any Green's operator can be fully broken down to basic operators if one can 
factor the differential operator T into first-order factors. Having a fundamental sys- 
tem ti, . . . , tm+n = Si, . . . , Sjn, . ■ . , T^n with succcssive Wronskians Wk{u) ^ {k = 
1, . . . , m -1- n) for u £ R+, such a factorization of T can always be achieved by well- 
known techniques described for example in Eqn. (18) of [15]; see also [17] and [23]. 
Using this factorization, we can break down G in a way similar to (3.3) except that 
the must be replaced by more complicated operators based on A and Si, similarly 
the by suitable operators involving B and rj. We assume m,n > throughout 
for avoiding degenerate cases. 

Proposition 3.2. The Green's operator of (3.1) is given by G = GsG,. where 
Gs — As-^ ■ ■ ■ Ag^ and G^ ~ (— l)"i?ri ■ • ■ -B^^ with 

^ As- ^ A^^ forl<i<m, 



Bt^Br =^ii2_5i£i^ form + l<i<m^ 

setting = 1 for convenience. 



Proof. We employ the factorization T ~ Tr^ ■ ■ ■ Tr^^Ts,^ ■ ■ ■ Ts^ , with the first-order 
operators given by 

Tt.=:^^D-^ forl<i<m, 

Tt.=Tr. ^ '}^Ll^D-^ for 771 + 1 < 7 < 771 + 71. 

It is then clear that G = Ag^ ■ ■ ■ Ag^^ (— ) • • • {—Br^) is a right inverse of T since 
both A and —B are right inverses of D. It remains to show that $ = Gg satisfies the 
boundary conditions. Differentiating $ fewer than m times results in an expression 
whose summands all have the form h ■ {A - ■ ■ g) ioi some functions h; evaluating any 
such summand yields h{0) ■ ( Jg ■ • • g) = 0, so the homogeneous initial conditions are 
indeed satisfied. For showing that the stability condition $(oo) = is also fulfilled 
we write $ = Ag-^g with g = Ag^ ■ ■ ■ Ag^^Grg. Then $ = siAs^^g and hence 

$(oo) = si(oo) j'^ si{u)~^g{u) du = 

because si(cx)) = and the integral is assumed to converge. □ 



Note that we assume, in the above proof and henceforth, that all forcing functions 
are chosen so that all infinite integrals have a finite value (this will be the case in all 
the examples treated here). This is also the reason why the rj are incorporated in 
B operators rather than in A operators as for the s^. Since we want to focus on the 
symbolic aspects here, we shall not elaborate these points further. 

Spelled out in detail, we can now write the Green's operator of (3.1) in the 
factored form 

G = ci ^ C2 ^ C3 • • ■ c„+„-i c„+„ (3.4) 

where Ci is A for 1 < z < 777 and —B for 777 + l<7<777-|-77. Although this brings us 
already some way towards a closed form for <I>(7i), we would like to collapse the m + n 
integrals of (3.4) into a single integration, just as we did in [2]. 
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To start with, assume for a moment that we did not have any unstable sohitions 
so that the fundamental system is only si, . . . , s„i. In that case we must dispense with 
the stability condition, imposing only the homogeneous initial conditions in (3.1). The 
Green's operator consists only of A operators, without any occurrence of B. In this 
simplified case, how can one collapse the m integral operators Ci, . . . , Cm = Am (3.4) 
by a linear combination of single integrators (multiplication operators combined with 
a single A)7 The answer is given by the usual variation-of-constants formula, which 
can be rewritten in our operator notation as follows [16]. 

Proposition 3.3. // si, . . . , s„i is a fundamental system for the homogeneous 
equation T<I> = 0, the Green's operator of (3.1) is given by 

Gs = s^A^ + --- + s,„A^, (3.5) 

where Wm is the Wronskian determinant o/si,--- , Sm and dm,i results from Wm by 
replacing the i-th column by the m-th unit vector. 



In other words, $ = Gg is a particular solution of T$ = g, made unique by 
imposing the initial conditions $(0) = $'(0) = . • . = $(™~i)(0) ~ 0. In our case, the 
stability condition $(cxd) = follows because 5^(00) = for all i = 1, . . . , m. But note 
that (3.5) is valid for any fundamental system si, . . . , of T, yielding a particular 
solution for the initial value problem (meaning (3.1) without the stability condition). 

Let us now turn to the general case, where the fundamental system ti, . . . , t„i+n 
consists of TO > 1 stable solutions si, . . . , Sm and n > 1 unstable solutions ri, . . . , r„. 
In that case the Green's operator has a representation analogous to (3.5) except that 
we need B operators in addition to A operators and we have to include definite inte- 
grals F for "balancing" the B against the A operators. 



Theorem 3.4. Define the constants 

CtiJ = rfj,m+j(0)/Wm+j-l(0) (3.6) 

for j = 1, . . . , n and i ~ 1, . . . , to + ri; the functions Oj ~ ai,j si + • ■ • + am.j Sm 
for j — I, . . . ,n; and the functions ai, . . . , a„ by the recursion ai = Oi, dj = aj — 
Oim+i.j di — ■ ■ ■ — am+j-i,j flj-i- Then the Green's operator of (3.1) is given by 

G-E^^^^^"E«.^^^^^' (3-7) 
where Ci is A for 1 < i < m and — _B for m + 1 < i < m + n. 



The proof of this result is given in Appendix A. There is a more explicit way of 
specifying the sequence of functions di, . . . ,d„ occurring in Theorem 3.4. 

Proposition 3.5. The functions dj in Theorem 3.4 can be computed by solving 
the system Td = a, where T is the lower triangular matrix with entries 

am+k{j) for j > k, 
Tjk ^ {I for j ^ k, 

otherwise, 



8 ALBRECHER, CONSTANTINESCU, PALMOWSKL REGENSBURGER, ROSENKRANZ 



while a and a are respectively columns with entries Oi, . . . a„ and Oi, . . . , a„ . Hence we 
have explicitly aj = dctTj/detT, where Tj is the matrix resulting from T by replacing 
its j-th column by a. 

Proof. We have am+i.j ai + ■ ■ ■ + am+j-ij 5,j-i + hj = aj, for j > 1. by the 
definition of the aj. But this is clearly the j-th row of the matrix Ta, while the 
recursion base oi = ai provides the first row. The explicit formula is an application 
of Cramer's rule. □ 



In either form, the functions ai, . . . , a„ can be readily computed from the given 
fundamental system si, . . . , Sm, ^i, . . . , r„, and the representation (3.7) provides a 
closed form for the Green's operator of (3.1). 

4. Asymptotic results for the renewal risk model. In the sequel, we will 
write fc(u) ~ l{u) if lim„^oo = 1 for some functions k and I. Assume that 
both the interclaim distribution and the claim size distribution have rational Laplace 
transform, i.e. their densities satisfy the ODE (2.1) and (2.3), respectively. Assume 
that the solutions of Equation (2.4) are of the form ti{u) ^ 

ti{u) ^ exp {Ar]i{u)} , i = l,...,n + m, (4.1) 

with 

rit{u) yi + i = l,...,n + m (4.2) 

u 

and 

y„i < ■■■ <yi <0 < yrn+i < ■■■ < ym+n (4.3) 
(so the rji are not asymptotically equivalent). 

Remark 4.1. Note that for the premium functions p{u) = c+eu, p{u) ~ c+ i^^^ 
and p{u) ~ cexpe/u, the corresponding ti fulfill the conditions (4.1)-(4.3). For 
m = n = 1, a more detailed analysis is presented in Section 5.3. 

Define the constants hk = — Y^^=i '^jk F '^^^+3.^+^ wifh appearing in (3.2) 
and ajk as defined in (3.6). For a permutation ip on {1, . . . , m + n} we define 

where sgntp denotes the parity of ip. 

Theorem 4.1. If g{u) ^ e""" for v > -j/i, then under (4.1)-(4.3) the asymp- 
totic expansion 

m+l 

$(u) « ^ Mu) (4.4) 

i=l 

holds, with 'i?i(u) = hiSi{u) (i = 1, . . . , m) and 

ra+n 

^m+i{u) - V — '—g{u). 

V i + V 
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This is equivalent to saying that hni„_j.oo — 1: fc = 1, . . . , m. 

Prool Note that by (4.1)-(4.2), tf' (u) ~ yfe^'''^"^ Using (3.7) and the Leibniz 
formula for the determinant, after some calculations one gets that expansion (4.4) 
with '&k{u) = lktk{u) and 



■dm+l{u) = ^ ■Kiti{u)C^ — {u). 
■ 1 

Using THopital's rule completes the proof. □ 

5. Compound Poisson risk process with exponential claims. Let us now 

focus on the case of a compound Poisson model (exponential intcrclaim times with 
mean 1/A) with exponential claim sizes with mean fi and a generic premium function 
p{u). The differential equation in (2.4) has order two in this case, so we expect to 
have one stable solution s and one unstable solution r. In fact, here we can relax the 
notion of an unstable solution, allowing any function where 

r(cxj) — lim r{u) 

exists and is different from zero (so the limit does not necessarily have to be infinity). 
The reason for this extension is that the basic argument for the ansatz 

carries over: Every solution of (2.4) must be of the form (5) since r{u),s{u) forms a 
fundamental system. But then the stability condition $(oo) = can only be satisfied 
if 7r = because we require s{oo) = 0. This is why the form (3.2) is still justified in 
the special case n = 1 with 71 — ^s- But note that this argument fails when there 
are more than two unstable solutions since they can cancel out unless we take some 
further precautions (e.g. requiring them to be of the same sign). 

5.1. Closed-form solutions for generic premium. For a discount factor 5 > 
0, the expected discounted penalty functions satisfies the second order LODE 

{D + fi) {-p(u)D + ^ + A) $(u) - A^$(u) = \{D + /i) uj{u). 

Expanding the operators, the equation is equivalent with 

{~p{u)D^ ~ {np{u) + p'{u) - X-d)D + d» $(u) = X{D + fi) uj{u). 

Assuming that p{u) ^ for all u > 0, this is further equivalent to 

+ + ^ - ^) - 4^) ^u) . ,(.), (5.1) 

with g{u) = ~j^iD + fj,) uj{u). Furthermore, we assume that p{u) is chosen in such 
a way that the associated homogeneous solution has a fundamental system s, r with 
one stable solution s and one unstable solution r with Wronskian w = IU2 = sr' — s'r 
nonzero on . Then the Green's operator for the boundary problem for the Gcrbcr- 
Shiu function $ is given by Theorem 3.4 with si = s and ri ~ r, namely 

G9{u) = ( - sin) J^^- riu) ^ + / ^S)-^^-) (^.2) 
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For calculating the full expression 

=7s(m) + G5(w) (5.3) 

we have to determine the constant 7. Evaluating the integro-differential equation 
(2.2) at zero, one obtains 

-c$'(0) + {X + 5) $(0) = Aw(0) 

and therefore 

Ac.(0)+c(Gg)'(0) A^O)+c '-(°>'%;'(°)-(°) /rgag(.)rf. ,,4, 
^ (A + 5)s(0) -cs'(O) (A + (5)s(0)-cs'(0) 

for the required constant. For 5 — 0, the LODE (5.1) is of first order in and 
its associated homogeneous equation has an unstable solution r{u) ~ 1 and a stable 
solution 



(5.5) 



(cf. [4]). For the fundamental system s,r, the Wronskian is just w = W2 = — s', and 
the Green's operator (5.2) specializes to 



00 



Ggiu) = (.(.) M~W>1 9i-)dv 

while the constant in $(u) = 7 s{u) + Gg{u) is now given by 

A^(o)-p(o)^ri^g(.)rf. 

^ As(0) -p(0)s'(0) 

Thus the Gerber-Shiu function can be written generically as 

Ac^(O) - p(0)^ 4t\ g{v) dv 
^' Xs{0)-p{0)s'{0) 



in terms of s(u). 



Remark 5.1. For S = and w = 1, one has g = and i/'(u) = 7s(w), recovering 
(1.3) for the ruin probability. 

5.2. Closed-form solutions for some particular premium structures. A) 
Linear premium: As discussed in Section 1, the linear function p{u) — c + eu can be 
interpreted as describing investments of the surplus into a bond with a fixed interest 
rate e > 0; see for example [21]. For S > and p{u) = c + eu, we can compute a 
fundamental system for the second-order LODE 

(D'^+(n+^- ^) D - <P(u) = --^ (D + a) uiu) 

\ ' V ' c+eu c-\-eu } c+eu } ^ ' c+eu ^ ' / V / 
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in the form 

s{u) ^U{l + 1,^ + 1, fiu + f ) {eu + c)^ cxp{~^^u), 

\+s ^^-^^ 
r{u) = A'/(f + 1, ^ + 1, (111+ if) {eu + c)~ exp(-^u), 

where M and U denote the usual Kumnicr functions as in [1, §13.1]. For u —¥ oo, the 
estimate in §13.1.8 yields 

s{u) = {eu + cf'-^ cxp(-Aiu) (l + ^ 0, (5.7) 

while the estimate in §13.1.4 yields 

r(u) = K2 {eu + c)^'' + ^ K-2 [eu + c)^/' +0{{eu + cY'^--^) ^ oo, (5.8) 

where Ki and K2 are some constants. Hence s is indeed a stable and r an unstable 
solution. Using §13.1.4 one derives the Wronskian 

-2 = '-^'-^{it^''"{ue + c)(^+*-)/-i exp(-M- + f ). 

Substituting these expressions in (5.2), we end up with 



Ggiu) = raZtniW e iff''''^' + -)''^''^' -M-i^u f )x 

( - Uiu) / Miv) - M{u) / U{v) + ^ U{u) / Uiv)) g{v) dv, 



where U{u) and M{u) are Kummer functions appearing on the right hand side of 
(5.6). This jointly with (5.4) is sufficient to determine the discounted penalty function 
in (5.3). 

B) Exponential premium: In general, an exponential premium function leads 
to intractable results. However, for 

p(m) = c(l + e-") 

the probability of ruin can be worked out from the expression in Section 5.1: 

^(^) ^ _ (1 + m; 1 + ^; e« + l)(ie" + i)^ 



2^^^(1 + ^,1 + 1; 2 + I; 2) 

where F{a, b; c; z) = 2-^1(0, b; c; z) stands for the hypergeometric function, see e.g. [1]. 
C) Rational premium: For a basic rational premium like 

P{u) = c + 

1 + u 

the exact symbolic form for the probability of ruin can be computed up to quadratures, 
namely 

A(C + 1) V';' g-n(cp-A)/c(^ _^ ^ ^y{X+c^)/c^ ^ ^j)^^ 



1 + A(c + l)Vc^ e-"(=A'-A)/c(c + 01* + l)-(A+c2)/c^ (1 + u)dv 



D) Quadratic premium: For the quadratic function p{u) = c + u^, the 
probability of ruin can be determined as 



^(u) = — - 

l + A/o°°e-(-^'^'-'=t^"(^/v^)+^^\^)/\^/(c + x2)da;' 
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5.3. Asymptotic results for generic premium. Assume the LODE 

+ Cn-iiu) + • • ■ + co{u) $ = (5.9) 

has complex coefficients Ci{u) contimious on and define its characteristic equation 
as 

+ Cn-iiu) + ■ . • + caiu) = 0. (5.10) 

Then the asymptotic behavior of solutions of (5.9) for u cx) essentially depends on 
the behavior of the roots j/i (u) ,...,?/„ (w) of (5.10) as u oo (see e.g. [7, §5.3.1, 
p. 250]), which will be exploited below. 

5.3.1. Probability of ruin. When 6 = and w = 1, the expected discounted 
penalty is the probability of ruin. For this quantity we have the following asymptotic 
estimate (we use the convention p{oo) = lim„^oo p(w)): 

Proposition 5.1. 
1. If p{oo) = c, where c is constant, then 



u — > oo. 



u oo. 



^(«)^^7exp(-^i^i + A^''^ 
2. If p{oo) = oo, then 

Proof. Integration by parts in (5.5) gives 

o^p{-t^v + xJ^ ^)rf.--exp(-M. + Ay^ ^). 
with s(0) = jh{^) — J, s'(0) = — kyr, where h denotes the Laplace transform of 



p(0)' 

h{u) = exp(A/J' ^). Letting f{u) = i exp (-^iv + XJ^ ^) dv, one gets 



^{u) = 7s(u) = t,f{u) ~ fiu) = C*A^)f{u), 

du 

^i*'^ ^ = As(o)-p(o)s'(o) - th^t = (l + ^) • To prove the first part 

of the proposition we need to show that V'(^) = ^l^i^f'{'^){^ + "(l))- According 
to our previous observation, ^•^J|^^7/|„|"'' = i+p(^u) ' '^^^ich completes the proof using 
I'Hopital rule. The second part can be proved similarly. □ 

5.3.2. Expected discounted penalty. We consider two cases of premium func- 
tions: 

PI. the premium function behaves like a constant at infinity 

p(oo)=c, p'{u)=o(^^- (5.11) 

or 
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P2. the premium function explodes at infinity, ^(oo) = cxd as 



= c + ^£iu'', oO. (5.12) 

i=l 

The first case is satisfied by the rational and exponential premium functions. The sec- 
ond case is satisfied by the linear and quadratic premium functions (see Section 5.2). 
Consider first the homogenous equation, (5.1) with g = 0, i.e. equation (2.4) with T 
given in (2.5) with 

s^l p'{u) x + s 

co{u) = ci(u =//+— — 

p{u) p{u) p{u) 

After tedious calculations one can check that in the case (5.11) we have 

r ci(u) 
lim < oo 

Vco(u) 

so that Conditions 1) and 2) of [7, p. 252] are satisfied. In the second case (5.12) we 
have 

lim ^=0 

and Conditions 1) and 2') of [7, p. 254] hold. From [7, p. 252] we hence know that 
for both cases (5.11) and (5.12) the asymptotic behavior of the solution of (5.1) is 



U{u)r^e^p\l {g^it) + g\'\t))dt\ , ^ = l,2, (5.13) 



where 



Qi = 



+ Pin) Pin) ) y[t^+ Pin) ' ' ' + ^ ' 



4 

) p(«) / pM 



and 



92 = ^ 

are the negative and positive solution, respectively, of the characteristic equation 

/ p'(u) X + S\ 5u , , 

x'+[fi+i^-—-)x-^^0. (5.14) 
V p{u) p(u) J p[u) 

Here and 0^2^^ arc defined by 

sf >(„) 
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Remark 5.2. Note that if the premium function p{u) satisfies conditions (5.11) 
and (5.12), the solutions ti{u) will be of the asymptotic form (4.1), where rji = Qi + Q^^'^ ■ 
In order to complete the asymptotic analysis of $(m) for large u, recall that the Gerber- 
Shiu function $ is given by ^{u) — ^s{u) + Gg{u), for a normalizing constant 7. 

Theorem 5.2. Under the assumptions (5.11) and (5.12) regarding the premium 
function, the asymptotics of the Gerber-Shiu function are described by 

$(u) ~ his{u)+Kig{u), 

with the exception 

$(u) ~ his{u) + K2ug{u) 

for the case (5.12) with 1 = 1. Here /ii = 7 — s'^i]) 9 {''-') dv. 

Remark 5.3. Moreover, for the particular examples considered here, the struc- 
ture of s (and r) is indeed that of the form (4.1)-(4.3) that we had to impose as a 
condition in the more general framework. 

Proof. First note that for (5 = 0, Gg{u) = 0, and thus $(w) has the same behavior 
as the probability of ruin V'(it) = ^s{u). Evaluating the expression (5.13) at J = 0, 

s{u) - e"^"+^ ^» ^ (Mu) + P'iu) - A)"^ 

leads to the classical result regarding the probability of ruin (1.3). For 6^0, one 
needs the asymptotic behavior of Gg{u), which based on (5.2) can be reduced to 
analyzing the asymptotic behavior of 



q{u) = -s{u) 1^ ^g{v)dv~r{u) j ^ g{v)dv, 

s{v) ^ 
s' (v) - 



since the term s{u) p7^5(^) ^'^ behaves as s{u) at infinity. After rewriting 



, , IoT^9iv)dv J^Il^g{v)dv 
= ^ 4 

s{u) r{u) 

and expanding the Wronskian, one can apply I'Hopital rule and see that as u — > 00 
(after some algebra) 

r {u) s [u] ^ ' s [u] r [u) ^ ^ ' 

q(u) ~ ^"5^ -Ti^-TM _ (5^15) 

s{u) r(u) 

Using Fedoryuk's asymptotic expressions (5.13) one more time, one can perform the 
analysis along the two cases introduced here. It easy to check that in the first case, 
PI, we have 

s(u) ~ exp {— fciu} , r(u) ^ exp {— fc2ii} , (5.16) 
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where ki = '-^ and fc2 = . Ihus, 

as lim„^oo 7^ = fci and lim„^oo 7^ = then 

h{u)^Kig{u), u — !■ cx), (5-17) 
where Ki = , ,^)t^'^ } \ = , ^ The second case, P2, is more complex, 



producing more intriguing asymptotics. One can show that in this case 

s{u) - hi u^e-'™, 

with /3 G M. Note that for I = 1, ei = e one recovers the asymptotics (5.7). One can 
also check that 

lim — — = -/i, 

u^oo S\U) 

whereas 

r'iu) 1 

— ~ - for I = 1, and r(u) ~ 1 for / > 1, 
r(u) u 

producing respectively 

q{u) ^ ug{u) and q{u) ~ (5.18) 

□ 

Example 5.3. Consider a compound Poisson risk model with premium functions 
described by assumptions (PI) or (P2). Let the penalty be a function of the surplus 
only, w(x,y) ~ e^'^^ . Since we are in the exponential claims scenario, 

/•OO 

g{u) ^ Xfi(D + n) / w(ii-y)e-^^d2; = Ai/e-(''+^'". 

J u 

Thus, for a linear premium function, 

with /3 ~ X/e — 1, whereas for all the other premium funtions in the class considered 
here. 



$(u) - /iiM^e"^" + Ai/e" 



u 00, 



with /? G K. 

Example 5.4. When p{u) = c expe/u, one has a differential equation with almost 
constant coefficients, 

Z?^ + (yu - 4j - exp -e/u)D - — exp -e/w ) <J>(u) 

c c J 

= -- exp -e/u{D + fi)uj{u). (5.19) 
c 

This is an equation of form (5.9), with coefficients satisfying 

Ckiu) = ak + okiu), k^l,2 (5.20) 
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with Q!fc constant and \ak{u)\du < oo. Here ai{u) ~ + ^^^a{u) and ao(w) = 

^a(u), with a(u) ~ X]fc!Li ^ ^1^] thus \a{u)\du < oo, and similarly for ag and 
ai. From [5, Th. 8.1, p. 92] (see also Problem 32, p. 105 there) we can hence conclude 
that the homogeneous equation has a fundamental system with asymptotics 

s(w) = e'""(l + o(l)) and r(M) = e"" (1 + o(l)) , 

where cr and p are solutions of the equation 

X + 6\ 5ii 

]x = 0, 

c / c 

with Re(o-) < and Re(p) > 0. Note that these solutions coincide with the one of 
the constant premium case. Consequently, one has the same asymptotic behavior as 
when the premium rate is constant. 

6. Conclusion. We have provided a symbolic method and a conceptual frame- 
work for studying boundary value problems with variable coefficients as they appear 
in modelling the surplus level in a portfolio of insurance contracts in classical risk 
theory. The approach presented allows a detailed analysis of the asymptotic behavior 
of the solutions of these equations under a set of conditions. For the specific case of 
the compound Poisson risk model, these conditions were made more explicit in terms 
of conditions on the form of p{u). Moreover, several new closed- form solutions were 
established within this framework. 

Appendix. Proof of Theorem 3.4. 

The proof of Theorem 3.4 hinges on the following technical lemma on Wronskian 
determinants. 

Lemma A.l. We have (^liti)' = -d^ for l< i < k < m + n. 

Proof. We have to show di^k+i w'i.—d[ j.^-^ Wk = di^k Wk+i- We note that all expres- 
sions in this formula are certain minors of the Wronskian matrix W for <i, . . . , tm+n- 
So let us write VTji' foi' the minor of W resulting from deleting the columns in- 
dexed ii, . . . ,ii and the rows indexed ji, . . . Then we have Wk = M^fc+i , dj.fc+i = 
{-iy+''+^Wl^^, diM = (-l)*+''W^fe;tti: with the derivatives 4 ^+1 = {-iy+''+^Wl 
and = W^~^^. For the latter, we use the fact that a Wronskian determinant can be 
differentiated if one replaces the last row by its derivative; see for example [10, p. 118] . 
Multiplying by (-1)'+'^', it remains to show WiW^+^ - Wi^^W^+^ = W'^IX\ ■ det W. 
But this is a classical determinant formula of Sylvester; see for example [20, p. 1571] 
or Eqn. (4.49") in [9]. □ 

The preceding lemma is the key tool for removing the nested integrals in (3.7). For 
seeing this, note that it can be read backwards as giving the integral oi di^k Wk+i/wJ,. 
In conjunction with certain operator identities taken from [18], this allows us to col- 
lapse expressions of the form A - ■ ■ A or B ■ ■ ■ B oi, at the interface of the two blocks, 
A---B. 

Proof. [Proof of Theorem 3.4] Note that the case n = reduces to Proposition 3.3, 
so we may assume n > in the sequel. We know from Proposition 3.2 that 




G^{~lTA,,---A,^Br,---Br., 
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using the notation employed there. Based on this factorization, we prove (3.7) by 
induction on n. In the base case n = 1, applying Proposition 3.3 again yields 



rn 



i=l 



drriA 1 B 



and Lemma A.l gives {dm+i^i/wm)' for the expression in parentheses. Now we em- 
ploy the identity AfB = A{J^f) + {j^f)B, for arbitrary functions /, from [18]. 
Substituting the expression in parentheses for /, this gives 

A + ^^^^ B-ai^F, 

SO we end up with 

G y f.s, A ^^^^ + B-^ -ai,s, F-^^) . 

j=l 

In the middle, we fector out si dui~\-i.i7 which equcils — vi drn-\-i,m-\-i 

as one sees by 

replacing the last row in Wm+i by the first and then expanding along that last row. 
But dm+i,m+i = Wrm SO the middle sum simplifies to ri {—B)dm+i,m+i/wm+i and 
may thus be incorporated into the first sum. In the third sum of the above expression, 
we factor out J2i "^i.i = ^i- Thus we obtain finally 

m+l 



which is the desired formula (3.7) for n — I. Now assume Equation (3.7) for n; we 
prove it for n + 1 . Using the induction hypothesis we obtain 

G = (-1)" An • • • A,,^Br, ■ ■ ■ i~Br„^, ) 



m+n 

' i=l j=l 

UCi{- di^m+n '"^2^""^') 



m+n 

E 



+ 

n 

As before, we see that Lemma A.l, with n + 1 in place of n, can be applied to the 
expressions in the two parentheses, yielding (dm+n+i.i/wm+n)' for the former and 
j/wm+n)' for the latter. In addition to the identity for AfB used for the 
base case, we need now also the related identities BfB = (J,^/) B — B (J^/) and 
FfB = F (J^/) also to be found in [18]. Whe n we substitute for /, these identities 
take on the form 
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for the first expression and 

FfB = F dm+n+l,m+j I Wrn+n — Ctn+l,m+j F 

for the second. Wc split the first sum above into the two sums 

y s, (A + 'h^i±^ B - a„+i , F) ^£2^, 

n 

E^./ ^^1 + ^ + 1. ni+j ^ _ _g '^m + 7i + l,m + j \ + 

In the lower-range sum 

m m 
i=l 

m 



2 = 1 

we can apply the same determinant expansion as before to obtain 

m n 

i=i j=i 

— fln+l ^ ; 



- ti Oi ^^^^^^^ + tm+n+l Om+„+i 

1=1 

n 

- y r, B ^^^^ - a„+i F ^^^^ ; 

in the upper-range sum we get 

n n 

E dm + Ti + l,m+j n "Jm + n , J. dm + ,i + l,m + j 

Combining the lower-range with the upper-range sum, the first sum within the latter 
cancels with the second sum within the former, yielding 

m+n+l 

Etj Cj "1+71+1,1 — ari-i-1 F — _ 
' ' ■U'm + n + 1 "^-^ Wm + n+1 

1 = 1 

Now let us tackle the second sum in the above expression for G, namely 



n 



+j_\ ;«„i+„ 
/ ii;„i+„+i 
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where the expression in parentheses is a„_|_i — a„+i by the definition of the Oj. Alto- 
gether we obtain now 



m+n+1 n+1 



which is indeed (3.7) with rt + 1 in place of rt. □ 



In concluding this Appendix, let us also mention that Theorem 3.4 is also valid 
if B is taken to be the operator J** with finite 6 G M rather than b = oo. The reason 
is that the operator identities from [18] are also valid in this case (and were actually 
set up for this case in the first place). 
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